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ABSTRACT 


Satellite soundings from the TIROS N Operational Vertical Sounder (TOYS) arc 
used to study the Experiment on Rapidly Intensifying Cyclones over the Atlantic 
(ERICA) Intensive Observation Period-2 (IOP-2), 13-14 Dee SS. TOYS data arc com¬ 
pared with dropwindsondes and coastal rawinsondes, and Spectral Grid Model (SGM) 
analyses of stability, geopotential height, and temperature. The impact of the first guess 
on the TOYS retrieval is studied by comparing four first-guess methods: 1) regression 
without surface data; 2) regression with surface data; 3) climatology with surface data; 
and 4) SGM analyses with surface data. Ship reports and moored and drifting buoys 
were objectively analyzed to obtain the surface data for the first guess. The retrievals 
with 6-h old SGM analyses as a first guess best captured the low level structure of the 
in situ soundings. The TOYS stability analyses defined the rapid cyclogcnesix environ¬ 
ment, in advance of tire developing cyclone, l'hc TOYS height and temperature analyses 
successfully described the structure of the developing cyclone, in reasonable agreement 
with experimental data and synoptic-scale analyses. 
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I. INTRODUCTION 


Historically, accurate weather forecasts have been limited by the availability and 
accuracy of observations. Even with increases in computer speed and the finer resol¬ 
ution of current numerical models, forecasts are still limited by the amount and quality 
of observed data. Over large data-sparsc areas, such as oceanic regions, important at¬ 
mospheric properties, such as temperature, pressure, humidity, and wind velocity, are 
not adequately measured. These limited data come front coastal rawinsondes, ship and 
buoy reports, and aircraft observations. An additional concern is that the middle part 
of the troposphere (about 5(.>0 mb) is infrequently sampled by these data, although con¬ 
sidered by meteorologists to be a key level for weather forecasting. 

One observation source that is capable of retrieving global temperature and 
moisture profiles is satellite soundings. Satellite-derived measurements of temperature, 
gcopotential height, and ozone have been available since the 1970 s. Their use in ini¬ 
tializing numerical forecasts has been steadily increasing. Over the data-sparsc southern 
hemisphere, satellite soundings are a primary observation source for NMC's current 
global forecast models (Dcy 19S9). 

Satellite soundings are important for military applications as well. In wartime when 
some data sources may not be available, satellite soundings can give forecasters a dy¬ 
namical perspective of the atmosphere over a much larger region than would be avail¬ 
able from a partial rawinsondc network. Particularly over the oceans, satellite sounding 
data could be especially valuable for military operations. The L’.S. Navy plans to im¬ 
plement real-time satellite sounding retrievals with the third phase of their Tactical En¬ 
vironmental Support System (TESS-3) in the 199U's. They will collect data from the 
NOAA polar-orbiting satellites to use as a nowcasting aid at sea and coastal locations. 
These data will include profiles of temperature and moisture from the Tiros-N Opera¬ 
tional Vertical Sounder (TOYS) onboard the NOAA satellites. The software available 
for retrieving TOYS soundings, during the early development phase of TLSS-3, was used 
in this thesis research. 

The purpose of this thesis is to evaluate the accuracy and synoptic utility of TOYS 
soundings over a coastal and oceanic environment during a rapid cyclogcnesis event. 
Data from the Experiment on Rapidly Intensifying Cyclones over the Atlantic (ERICA), 
conducted from December 1 OHS to February 19S9, wili be used for the evaluation. Rc- 
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search will concentrate on Intensive Observation Period-2 fIOP-2). 13-14 December 
19SS. 

In this study, TOYS soundings for IOP-2 will be compared with coastal rawinsondcs 
and aircraft dropwindsondcs. TOYS geopotential height and temperature plots and 
TOYS potential temperature cross sections will be studied in significant regions of the 
ERICA cyclones. Of particular interest is the determination if TOYS soundings can 
describe regions of low static stability observed by dropwindsondcs and coastal 
rawinsondcs available during ERICA. Emphasis will be placed on synoptic evaluation 
of the TOYS soundings and associated horizontal fields rather than computing error 
statistics (such as rms differences). 

The TOYS-derivcd soundings were prepared from the Naval Oceanographic and 
Atmospheric Research Laboratory (NOARL) version of the Simultaneous Physical Re¬ 
trieval TOYS Export Package, version 4.0. March PASS, developed by the Cooperative 
Institute for Meteorological Satellite Studies (CIMSS) at the University of Wisconsin 
(Smith et al. 19S5). All software was run on a microYAX III in the Interactive Digital 
Environmental Analysis (IDEA) Laboratory at the Naval Postgraduate School (NTS). 

The thesis is organized in the following manner: Chapter 2 discusses different 
sounding methods and their advantages and disadvantages; Chapter 3 describes the re¬ 
sults of various TOYS retrieval methods; Chapter 4 gives a background on ERICA 
IOP-2 and a more detailed description of the software used to produce the TOYS 
soundings for this thesis; and Chapter 5 describes the results of this study. All figures 
discussed in the thesis appear in Appendix A. 



II. SOUNDING METHODS 


A. THE TOYS PACKAGE 

The TOYS consists of three sensors that measure the emitted radiance at the top 
of the atmosphere. They arc: the second version of the High-resolution Infrared Radi¬ 
ation Sounder (IHRS), the Microwave Sounding Unit (MSU), and the Stratospheric 
Sounding Unit (SSU). The HIRS contains 20 (19 infrared and one visible) channels to 
determine temperature, water vapor concentration, cloud height and amount, and o/one 
concentration in the troposphere and lower stratosphere. The MSU contains four 
channels in the cloud-penetrating microwave spectrum to supplement the HIRS over 
cloudy regions. The SSU contains three channels to determine the temperature profile 
of the upper stratosphere. Fig. 1 (NO AA 1981) gives the physical characteristics of each 
of the wavelength bands of the 27 TOYS channels. Each channel is sensitive to a dif¬ 
ferent wavelength and is dependent on the absorption characteristics of 
C() : . // ; 0. 0 : , CU and A\() in the visible, infrared, and microwave spectral regions. Ra¬ 
diance measurements are made by the three sensors and converted into temperature, 
ozone, and moisture profiles by various retrieval methods. 

For each channel, the level of peak energy contribution in the atmosphere comes 
from the vertical layers indicated by the peak of the weighting function, or the change 
of atmospheric transmittance with height. The weighting function weights the 
blackbody radiance from the atmosphere, so it is sensitive to changes in atmospheric 
temperature for a particular layer. The TOYS weighting functions lor all 27 channels, 
giving the layer of peak energy contribution, are shown in Fig. 2 (NOAA 19M ). 

IF STATISTICAL RETRIEYALS 

Yarious methods have been devised to retrieve the TOYS radiance measurements 
and convert them into soundings of temperature, moisture, and o/one. Smith et al. 
(1979 ) describe a statistical approach. A database of climatological temperature profiles 
(training data set) based on latitude and surface type (land or ocean) is used as a first 
guess for the final retrieved temperature and moisture soundings. The training data set 
is needed to derive regression coefficients for converting the observed radiances to tem¬ 
perature and water vapor. It also gives the retrieval method a physically reasonable 
temperature profile that could have produced the observed radiance values. The re¬ 
gression coefficients arc updated weekly 



The horizontal resolution of temperature and moisture soundings from tire statistical 
retrieval method is determined from interpolated 111RS and MSL scan spots, or fields 
of view (FOV's). Each TOY is about 25 km in diameter. Radiance values are limb 
corrected so all FOV's can be treated as if they were a* nadir. A 9 X 7 TOY array is used 
to prepare one TON'S sounding with a horizontal resolution of about 250 km. 

The statistical method requires such a large TOY array because it requires dear- 
column radiances to determine the most accurate temperature and moisture profiles. 
TOYS soundings are based mostly on infrared wavelengths since 20 of the 24 sounding 
channels arc from the IIIRS. Therefore, the soundings are very sensitive to clouds. 
Channels that saturate when the atmosphere is too cloudy are avoided in the final tem¬ 
perature retrieval because they would produce temperatures that arc too cold. A 
sounding is considered to be clear if at least four TOY'S are clear. If fewer titan lour 
clear I OY s arc found, the retrieval tries to determine dear-column radiances from 
partly cloudy scan spots by determining the differences in cloud amount of adjacent 
pairs. If at least four adjacent pairs are found in which cloud differences can be dis¬ 
cerned. the sounding is judged to be partly cloudy. If fewer than four adjacent pairs are 
found, the sounding is determined to be cloudy, and the temperature retrieval is at¬ 
tempted with the MSL and stratospheric levels of the 111 RS only. 

The retrieval then produces vertical profiles of temperature, water vapor, mixing 
icio, gcopotcntial thickness, and total ozone concentration using regression methods. 
This statistical retrieval method was used operationally by NOAA from 1979 to 1989 for 
NOAA satellite soundings. 

C. PHYSICAL RETRIEVALS 

Another approach for temperature and moisture retrieval solves the radiative trans¬ 
fer equation (RTF) for brightness temperatures from the observed radiances. The so- 
called "physical retrieval method" described by Smith ct al. (1983) accounts for terrain 
elevation, surface emissivity. and surface temperature in solving the RTF. The physical 
retrieval uses a 3X3 array of interpolated I1IRS and MSL scan spots hr each sounding. 
As in the statistical retrieval, each FOY is checked for cloud contamination. Each 9- 
FOY box gives a sounding with a horizontal resolution of about 75 km. 

Lnlikc the statistical retrieval method, a statistical set of coincident rawinsondc 
profiles is not needed. Instead, first-guess temperature and moisture profiles must b: 
specified by one of three options: a climatological profile, a regression estimate from the 
observed radiances in the MSI.’ channels and stratospheric 111RS channels, or profiles 



produced by a numerical forecast model. The final temperature profile is determined by 
an iterative formula that requires first-guess temperature and brightness temperature 
profiles, weighting functions for 15 mandatory levels in tire troposphere and stratosphere 
interpolated from all the I11RS and MSU channels, and the calculated brightness tem¬ 
peratures from the RTE solution. 

The weighting functions required by the iterative formula are derived from calcu¬ 
lated transmittanees. The transmittanecs arc derived statistically using a set of re¬ 
gression coefficients that relate atmospheric absorption characteristics to temperature 
and moisture. The coefficients are generated for each of six equally spaced scan angles, 
alleviating the need to limb-correct the radiance measurements. They arc constant for 
each satellite, but can be changed if the sa eliite characteristics change (ie. sensor tem¬ 
perature change, etc.). 

The iterative formula is solved for temperature at 15 mandatory levels. Iteration 
stops when the difference between the brightness temperature and guess brightness 
temperature from the previous step is less than .05 K. Total ozone concentration is also 
calculated with an iterative formula, using the observed radiance of the ozone-sensitive 
HIRS channel (Cll 9) as a first guess. The water vapor profile is calculated using a 
mean relative humidity solution (described by Smith ct al. 19S3) from the observed ra¬ 
diances in the I {IRS water vapor channels, and the final temperature estimates from the 
iterative solution. 

D. SIMULTANEOUS RETRIEVAL 

A more recent version of the physical retrieval method, described by Smith ct al. 
(19S5), solves for surface skin temperature, temperature profiles, and water vapor pro¬ 
files simultaneously. The interdependence of the three parameters is taken into account 
and the importance of the first guess is reduced. The radiative transfer equation is solved 
in perturbation form, with respect to a mean atmospheric state. The solution for tem¬ 
perature and moisture perturbations requires a matrix inversion by least squares method. 
Weighting functions, as described in the physical retrieval method, serve as the basis 
functions for the matrix inversion. 

The simultaneous physical retrieval method is more computationally efficient than 
the iterative physical method. Instead of an iterative formula, a two-step method is 
employed for solving the perturbation form of the RTE. In the first step, the first-guess 
profiles of temperature, moisture, and skin temperature with their corresponding 
weighting functions, are used to solve for brightness temperatures at each pressure level. 



The same types of first-guess options arc available as in the iterative physical retrieval 
method. The regression estimate uses 111RS CH's 1-3 (stratospheric channels), MSI' 
CH's 2-4, and H1RS CH’s 11 and 12 (sensitive to water vapor), along with an estimate 
of the surface skin temperature from the warmest window channel. The RTE is then 
solved once more for surface skin temperature, and temperature and moisture profiles, 
using all channels determined to be free of clouds. The method for retrieving ozone 
concentration is the same as that described for the physical (iterative) method. 

The problem of determining cloud-free areas and retrieving profiles in cloudy regions 
is the most limiting factor in retrieving accurate profiles of temperature and moisture. 
In the simultaneous retrieval method, cloud contamination is determined by comparing 
the observed channel radiances with the first-guess skin temperature, the visible 
rellectance frem IIIRS CH 2<>. and the radiance values of the three MRS window 
channels (S, IS, and 19). Once the cloud contaminated channels arc found, cloud 
amount and height arc specified by using absorption characteristics of CO : in adjacent 
channels and the C0 2 slicing technique (Smith ct al. 1985). 

A more recent version of the simultaneous retrieval method, described by Huang 
and Smith (1986). improves upon the retrieval of temperature and water vapor profiles 
in cloudy regions. This later version includes the simultaneous retrieval of cloud amount 
and cloud top temperature in the inverse matrix solution of the perturbation form of the 
RTE. Although the later version was not available for this thesis, it will be available for 
incorporation into TESS-3. Preliminary results from the latest method are described in 
Chapter 3. 

E. COMPARISONS BETWEEN THE RETRIEVAL METHODS 

The simultaneous and iterative physical retrieval methods, described earlier in this 
chapter, were designed to be portable. They have been sent to many diUcrent users 
around the world. The latest version, the Simultaneous Physical Retrieval TOYS Export 
Package (Smith et al. 19S5), version 4.0. was sent to NOARE in March 1988. 

The advantages of physical retrievals over statistical methods are tlnn the RTE is 
solved in the retrieval and. thus, physical processes arc used throughout. Also, the data 
do not need to be limb corrected (a possible source of error). Visible channel reflectance 
(from MRS Cl I 20) is accounted for in the calculation of surface temperature. And 
finally, cloud contamination of the MRS channels is accounted for to allow maximum 
use of the MRS channels in a cloud}' atmosphere. 
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Disadvantages of physical retrievals are that the transmittance coefficients arc con¬ 
stant and assumed to be accurate. Poor quality soundings will result if the coefficients 
are not accurate and operational processing of the TOYS data may have to stop until 
the coefficients arc corrected. Statistical retrievals avoid the transmittance calculation 
problem because of the continuous updating of regression coefficients using radiosonde 
observations. 
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III. BACKGROUND 


This chapter describes several studies that evaluated the ability of TOYS soundings 
to describe significant synoptic and mesoscale features defined by in situ data. These 
investigations illustrate that TOYS soundings provide an excellent supplement to 
radiosonde data in analyzing a variety of meteorological phenomena. 

A. STATISTICAL RETRIEYAL STUDY 

Smith ct al. (1979) evaluated the statistical retrieval method by compositing TOVS 
soundings and radiosonde reports over N. America (Tig. 3). They found that TOYS 
soundings were within about 2°C of the radiosonde observations. The largest differences 
in the satellite and radiosonde observations were at the surface and tropopausc. where 
the vertical resolution of the TOYS is not fine enough to resolve inversions in the tem¬ 
perature profile. 

In the same paper, Smith et al. also evaluated the ability of the TOYS statistical 
retrieval method to describe the environment of an intense storm over Texas in March 
1979. They related TOYS 300 mb geopotential heights to a satellite image of the storm 
and plotted rawinsondc winds observed 9 h before and 3 h after the satellite pass on a 
TOYS 500 mb geopotential height analysis (not shown). The TOYS 500 mb and 300 
mb analyses showed a trough over Texas between the rawinsondc report times. This 
trough led to heavy snowfall on the Tcxas-Oklahoma border. The rawinsondc reports 
before and after the satellite pass did not clearly delineate the storm: however, the 300 
mb rawinsondc winds did show a wind shift between observation times, supporting the 
development of the trough. 

Smith et al. also evaluated operational TOYS 1000-500 mb thicknesses over the 
Northern Hemisphere. XMC 1000-500 mb thickness analyses were compared with and 
without TOYS data. The plots were strikingly similar (Fig. 4). The largest differences 
were associated with the more intense troughs over the eastern North Pacific Ocean and 
eastern Europe on the TOYS analysis. In summary, the study showed that TOYS data 
have a strong positive influence on global and local forecasting. 

B. PHYSICAL RETRIEVAL STUDY 

Smith ct al. (19S3) used the Physical Retrieval TOYS Export Package to explore the 
Alpine Experiment (AEPEX) data in Europe during March 19S2. They compared TOYS 
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gcopotential thickness analyses with thickness analyses produced from radiosonde data 
by the European Center for Medium-range Weather Forecasting (ECMWF). A re¬ 
gression first-guess with surface data was used to retrieve the TOYS temperature and 
geopotential height profiles. The TOYS thickness analyses were prepared using a Barnes 
objective analysis with a 75 km grid spacing, while the ECMWF analyses were based on 
optimal interpolation with a spatial resolution of 200 km. Smith et al. found that the 
TOYS 1000-500 mb thickness analysis resolved a higher amplitude ridge between En¬ 
gland and Spain, and a deeper trough over the Mediterranean Sea (Fig. 5) for one of the 
ALPEX cases. This difference was partly due to the coarser horizontal resolution of the 
ECMWF analysis. However, in the low stratosphere the situation was reversed. The 
300-100 mb TOYS thickness analysis (not shown) revealed a weaker horizontal gradient 
than the ECMWF analysis. The ECMWF analysis showed a stronger trough and ridge 
in these upper levels. The cause of the poorer 300-1O0 mb TOYS analysis was attributed 
to coarser TOYS vertical resolution in the upper levels. 

This study also looked at the impact of the first guess and surface temperature 
options in the physical retrieval by comparing standard deviations between TOYS and 
radiosonde profiles for four different cases. The cases were: 1.) regression first-guess 
with surface temperature, 2.) climatology first-guess with surface temperature. 3.) re¬ 
gression first-guess without surface temperature, and 4.) climatology first-guess without 
surface temperature. The greatest impact of using the surface data was below 700 mb. 
When surface data were included in the retrievals, the 1000-700 mb error was signif¬ 
icantly less (1-2'C for the regression first-guess option and 1-2.5X for the climatology 
first-guess option). Without surface data, the error increased to 2-3 c C from 1000-700 
mb for the regression first-guess option and 2.5-4.5°C for the climatology first-guess 
option. In conclusion, this physical retrieval study showed results similar to the statis¬ 
tical retrieval study. TOYS data gave physically reasonable estimates of atmospheric 
structure when compared with "ground truth" radiosonde data. 

C. SIMULTANEOUS RETRIEYAL STUDY 

Smith et al. (1985) applied the simultaneous retrieval method to the same ALPEX 
data set. Objective analyses of temperature and dewpoint, derived from radiosonde ob¬ 
servations at 1200 UTC 4 March and 0000 UTC 5 March 19S2, were compared with 
TOYS analyses at 850. 700, 500, and 300 mb. 

This study first compared two first-guess options of the simultaneous retries al 
method: regression and climatology. Smith et al. found that the two first-guess moth- 




ods gave similar results for all levels except 3<X> mb. The regression first-guess retrievals 
were slightly better at all levels, except 3u0 mb, when compared with the ECMWE ana¬ 
lyses. However, in all cases, the differences were generally within 3°C cf the ECMWF 
analyses. The largest differences were found over the radiosonde data-void areas. Also, 
both analyses revealed the rapid intensification of a trough over western Europe between 
1200 ETC 4 March and 0000 ETC 5 March 19S2. The first guess did not have a major 
impact on the temperature retrievals in this study. 

Next, the simultaneous retrieval method was compared with the physical (iterative) 
retrieval method. The temperature profiles derived by the simultaneous retrieval were 
more accurate than the physical retrieval method by .5 C C at nearly all levels up to 250 
mb (Eig. 6). The RMS difference for the simultaneous retrieval method with a regression 
first-guess was generally the same as that with a climatology first-guess from 1000-500 
mb. The difference between the two was significant only between 300-500 mb and at 
200 mb. showing an improvement of 0.5-1.(EC for the regression first-guess. Smith et 
al. stated that the coefficients used for the regression first-guess helped the retrieval 
method derive a more accurate tropopause structure. 

However, the greatest difference found between the two methods was in the 
dewpoint results. Smith ct al. found a significant improvement in the dewpoint analyses 
using the simultaneous method. The simultaneous method also resolved larger and more 
meteorologically consistent dewpoint gradients. 

finally, the study explored the impact of surface data on the first guess in the si¬ 
multaneous retrieval method. The greatest effect of surface data on the retiieval was 
from 1000-700 mb (Eig. 7). Near the surface, the error decreased by half when surface 
data were used. In summary, this study showed that the simultaneous retrieval was an 
improvement over previous retrievals. It gave better results than the physical retrieval 
and was more efficient. The investigation emphasized that TOY'S soundings are useful 
and provide an excellent supplement to other data sources. 

D. SIMULTANEOUS RETRIEVAL WITH CLOUD PARAMETERS 

Huang and Smith (1986) revised the simultaneous retrieval method to include cloud 
amount and cloud top temperature in the matrix inversion solution of the RTE. In a 
2-stcp simultaneous method, the C0 2 slicing method (Smith ct al. 1985) is performed 
after the temperature and dewpoint retrieval to determine cloud parameters (height, 
amount, and temperature). 
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The study compared TOVS temperature retrievals with and without cloud parame¬ 
ters during ALPEX. The soundings improved during all comparisons of clear, moder¬ 
ately cloudy, and extremely cloudy conditions, but the greatest improvements were 
achieved in the extremely cloudy cases. RMS differences throughout the troposphere for 
the extremely cloudy cases (around 2°C to start with) decreased to an average of 1.3°C 
with the revised simultaneous retrieval method. 

E. OTHER RETRIEVAL STUDIES 

1. Goddard Laboratory for Atmospheres (GLA) physical retrieval 

Other groups besides the University of Wisconsin CIMSS have made progress 
using TOVS. Susskind et al. (1984) described a physical retrieval method developed at 
the GLA. The retrieval method is iterative and produces soundings with a resolution 
of about 125 km. The RMS statistics showed that the GLA physical retrievals were 
much better than the operational soundings for a first GARP Global Experiment 
(EGGE) data set. The largest errors were at the surface (2.8 K). decreasing up to 700 
mb (2.2 K). and remaining the same throughout the rest of the troposphere. 

The GLA retrieval also gives derived sca-surlace temperature fields, fractional 
cloud coverage, cloud top pressure, cloud top temperature, surface emissivity, and sea 
ice extent by interpreting the radiance observations of all 111RS and MSE' channels. 
Retrievals for all of these fields were found to be reasonably consistent with ground- 
truth data. 

2. Dependence of the first guess on the GLA physical retrieval 

Reuter et al. (1988) investigated the dependence of first-guess data on the GLA 
physically-based retrievals. Their study, like Smith et al., used the ALPEX data (1290 
L1C 4 March 1982 - 0000 ETC 5 March 1982). The retrieval used for the study modi¬ 
fied the Susskind et al. (19S4) physical retrieval by improving the moisture retrieval al¬ 
gorithm. In addition, experiments were conducted to improve the first guess if 
determined to be poor. Susskind and Pursch found that the thickness Helds showed a 
low first-guess dependence. This result agrees with the Smith ct al. (1985) study, 
However, the prceipitablc water fields did show a stronger first-guess dependence, espe¬ 
cially close to the surface. 

3. Using a forecast model-output first-guess to improve TOVS soundings 

Daniels ct al. (1989) described a study to improve the XLSD1S operational 
temperature and moisture retrievals by inputting NMC’s Medium Range forecast 
Model 6-h forecast as a first guess instead of the (then current) operational method of 
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using a data set of radiosonde temperature and moisture retrievals and satellite radiance 
measurements as a first guess. Their results showed a significant improvement over the 
operational method. 

4. Simultaneous retrieval over Australia 

Lc Marshall ct al. (19S9) studied the accuracy of the simultaneous retrieval 
method (Smith ct al. 19S5) over Australia. The retrievals used either a statistical 
retrieval-based or a model forecast first-guess for temperature and moisture profiles. 
The retrieval method has been operational since late 19S7. The study compared TO VS 
retrievals over the Australian region with NMC analyses. 

TOYS potential temperature cross sections. 1000-500 mb thickness analyses, 
ozone concentrations, temperature profiles, and temperature analyses were compared 
with in situ data. The largest RMS difference between TOYS and NMC temperature 
profiles occurred at the surface and tropopause (2.5 K). The study found a high corre¬ 
lation between TOYS and NMC 500 mb temperatures over Australia. TOYS potenti.il 
temperature cross sections and 1000-500 mb thicknesses were also very similar to the 
NMC analyses. Overall, the TOYS analyses compared extremely well with the NMC 
model data. 

The above studies have shown that atmospheric structure can be accurately 
described by TOYS data. TOYS structure supplements in situ data sources, such as 
radiosondes, by detection of mcsoscale features (both temporally and spatially) within 
the radiosonde network. TOYS data, using the Simultaneous Retrieval Method with a 
variety of first-guess options, and with and w ithout surface data, will now be applied to 
the I;RICA data set. 
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IV. ERICA: BACKGROUND AND DATA DESCRIPTION 


A. ERICA BACKGROUND 

As indicated in the previous chapter, the rapid cyclogencsis event investigated in this 
study is I OP-2 (13-14 December 19SS). Some preliminary, research has already been 
completed on this case. Chalf'ant (1989) showed that the early cyclogencsis period was 
characterized by multiple surface centers. The eventual rapid deepener formed as a wave 
along a pressure trough cast of Cape Hatteras at about 1SO0 ETC 13 December 1%S. 
It rapidly intensified over the next 12-24 hours, deepening more than 30 mb from 
0000-1200 UTC 14 December, as it tracked slowly eastward. Chalfunt found that 
NMC's Spectral Grid Model (SGM> was unable to accurately predict details of the de¬ 
velopment of this cyclone during IOP-2. However, longer range forecasts (45s-h) were 
successful in capturing the magnitude of the deepening. Miller (1990) showed that the 
Naval Operational Regional Analysis and Prediction System (XORAPS) forecasts for 
lOP-2. were quite sensitise to the initial analyzed fields. Miller's study also showed that 
NORA PS had difficulty resolving the multiple surface centers present in tire early lOP-2 
period. Additional data now available for further research include dropwindsondcs 
around the storm and TOYS soundings covering the whole ERICA region. 

The additional mcsoscale observations available during ERICA are proving critical 
in the diagnosis of rapid cyclogencsis. Unfortunately, those observations will not be 
available in a future operational environment. On the other hand, satellite soundings 
may be able to provide some of this key information over the data-sparse regions. If 
important features such as low static stability regions and short wave troughs and ridges 
can be accurately described by the TON'S data, this additional information can be in¬ 
corporated into analysis and forecast models to more successfully predict rapid 
cyclogencsis. 

B. DATA DESCRIPTION 

The dropwindsondcs and rawinsondes available during ERICA provided a useful 
data set to compare with the TON'S soundings. However, rawinsondes give a much 
better resolution than the TOYS soundings. The I ON'S only contains 12 111RS channels 
for determining temperature profiles in tire troposphere and lower stratosphere (up to 
100 mb). Recall that the TOYS weighting functions (Tig. 2t have a fairiy broad vertical 
width of approximately 300-400 mb. This means that the radiance in each ION'S 
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channel conics from a relatively thick layer of the atmosphere. Therefore, a point com¬ 
parison between a TOYS sounding and a rawinsondc may be misleading, especially over 
a rapidly evolving cyclogenctic environment. So, for this stud}’, additional comparisons 
were made between TOYS geopotential height and temperature fields. The ability of 
TOYS to physically represent the environment described by the SGM (not taking into 
account possible SGM errors) then can be assessed. 

As mentioned earlier, a cloudy atmosphere presents the biggest challenge in retriev¬ 
ing an accurate satellite sounding. The problem with retrieving TOYS soundings over 
cloudy areas is that the I HRS channels cannot sense through clouds. Because the TOYS 
contains only four MSU channels, vertical resolution in cloudy regions can be expected 
to be poor. However, the simultaneous method tries to account for this problem by 
using ail I HRS channels sensing above the cloudy later and determining a "reasonable" 
profile below the top of the cloud and adjusting it with the MSU-dctermincd profile. 
Unfortunately, if the top of the cloud is determined to be at the 700 mb level (MSL CI1 
2) or higher. TOYS cannot determine any structure in the critical low troposphere below 
700 mb. MSU Cl 1 I (window channel) can penetrate through clouds, hut will give an 
inaccurate profile in precipitating regions. 

C. NO A RL/.NTS TOYS RETRIEVALS 

The retrieval method used for this study was a version of the Simultaneous Physical 
Retrieval TOYS l:\port Package (Smith et al. 1985). version 4.0, restructured by 
Richardson and Haugen at NOARL. The NOARL version, redesigned in f ortran 77 
format, speeds up the retrieval and allows for easier upgrades. A description of the 
processing steps in the NOARL version follows. 

The retrieval process begins with a check of the warmest HIRS CI1 <N (IR II pin 
window channel) brightness temperatures (/.) in a 9-TOY box. Recall that each DIRS 
TOY is about 26 km in diameter, so the resultant resolution of an individual sounding 
(with MSU interpolated TOV's) represents about a 7 5 km volume. Next, the scan angle, 
local zenith angle, and surface elevation arc determined. Tor this thesis, surface topog¬ 
raphy with 111 km resolution was used. Higher resolution topography (18.5 km), 
available in the original export package, was not available in the NOARL version. This 
higher resolution topograph}’ is not critical for this study since most of the area of in¬ 
terest is over the ocean. Alter the elevation is determined, the retrieval then finds all 
CH 8 7fs within 3 K of the warmest TOY in the 9-TOY box. This determines which 
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FOV's to use throughout the retrieval process. It then averages all of the warm I OV 
temperatures for each 111RS tropospheric channel to represent one TOYS sounding. 

The next step of the retrieval method is the most critical of all. The retrieval requires 
first-guess temperature and moisture profiles and 1000 mb heights for the initial solution 
of the RTE. Three options arc available. First, it obtains climatological profiles of 
temperature and relative humidity (-40 levels in the troposphere and stratosphere). 1000 
mb heights, and a surface temperature extrapolated dry adiabatically from 1000 mb. 
The climatological values arc based on five latitude zones, summer and winter. The re¬ 
trieval interpolates between the zones and seasons depending on the date of the satellite 
pass. Using the hydrostatic approximation and equation of state, the retrieval also cal¬ 
culates a station pressure. The user can replace climatological profiles with model- 
generated temperature and moisture profiles, including surface temperature, dewpoint, 
and lOoO mb heights. Also, if a regression first-guess is chosen, the climatological tem¬ 
peratures are replaced with the regression estimates front 1IIRS CITs 1,2, and 3 
(stratospheric channels) and MSU (ill's 2, 3. and 4. 

As mentioned in Chapter 3, Smith ct al. (19X5) showed that the difference between 
climatology and regression first-guess options was not significant, but that surface tem¬ 
perature and dewpoints estimates were critical in the final retrieval. If surface data are 
available, the retrieval method accepts the estimate as "ground truth" and weights it 
higher in tire final retrieval. In that way. the TOYS radiance estimates of surface tem¬ 
perature do not contaminate good estimates of the low level atmospheric structure. 
User-input 1000 mb heights arc also accepted as ground-truth and are not adjusted in 
the final retrieval. The use of model-generated temperature, moisture, and 100o mb 
heights as a first guess was not discussed in Smith ct al.'s study but will be used in this 
thesis. In addition, the regression and climatology options will also be employed in this 
study. 

If no surface data are available, the retrieval tries to improve the initial estimates 
of surface temperature and dewpoint. 1 irst. it searches for the closest pressure level to 
the surface. If the station pressure is less than IOOO mb (high ground), it recalculates the 
surface temperature, dewpoint, and mixing ratio based on the maximum of IIIRS CH S 
and 10 (8.3 iiinll 2 0>). MSU CH 1 (microwave window), and the 1000 mb first-guess 
temperature. If a super-adiabatic lapse rate is found at the surface, it is reduced to 
dry-adiabatic between the surface and the closest pressure level. Finally, if the closest 
pressure level is not equal to 1000 mb. the retrieval extrapolates the temperature of the 
closest pressure level isothcrmally to 1000 mb. 
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The surface skin temperature is also required in the simultaneous retrieval method. 
It is derived from 1IIRS CM 8 or 18 (IR 4 nm window) or MSL Cl I 1. First, the re¬ 
trieval uses the solar zenith angle to calculate visible reflectance from I HRS Cl I 2 ( >, If 
the reflectance is less than 25%, then the Cl I IS 7j, is used as the skin temperature . 
Otherwise, the maximum between Cll 8 and MSL' CH 1 is used, to avoid the cloud-top 
temperatures seen by CM 8. 

The retrieval then attempts to find which IIIRS tropospheric channels aie too 
cloud-contaminated to use. It assumes a cloud level where the difference between the 
CM 10 T b (WO mb //,<?,) and the guess T b for each channel is less than 10 K. or wherever 
a temperature inversion is encountered. The matrix inversion is then performed to ob¬ 
tain temperature and moisture profiles. The weighting functions for I MRS CM 1 and 
MSU CIFs 2-4 serve as the basis functions for temperature and I11RS CM T and CM 
12 weighting functions serve as the basis functions for water vapor in the initial retrieval. 

After the matrix inversion, the skin temperature is updated and used in a repeated 
matrix-inversion solution (second step). Cloud height and amount arc redetermined and 
lilRS channels contaminated by clouds are rejected before the second step. However, 
the retrieval only rejects the shortwave channels (4.? um) or the longwave channels (15 
urn), but never both. The influence of 1I1RS channels below cloud level is reduced in the 
second step. Also, if the cloud height is determined to be above 450 mb. the final tem¬ 
perature and moisture soundings are output as "cloudy" and the second step is not per¬ 
formed. If the cloud height is below 43o mb. the second step is preformed to determine 
a new temperature and moisture profile. Cloud height and amount is again checked. 
If the cloud level is above 570 mb, the final retrieval is output as a "cloudy" sounding. 
Otherwise, the sounding is output as "clear". Geopotential heights for all mandatory 
levels are calculated using the hydrostatic approximation and equation of state. How¬ 
ever. if 1000 mb heights arc not available in the first guess, they should be used with 
caution, since a climatological 1000 mb height is used. 

For this thesis, the various first-guess options were studied to determine the best 
combination of choices. Four first-guess options were explored. First, the retrieval was 
run with a regression first-guess and no surface data to determine the accuracy of TOYS 
data in a region where no other data sources may be available, such as in a wartime 
environment. Next, the impact of surface temperature on the retrievals using both 
climatology and regression first-guess options was studied to determine how much im¬ 
provement. if any, could be made to the TOYS data. The surface data included addi¬ 
tional marine data (buoys and ship reports) available during FRK A. I he data were 
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objectively analyzed using the Barnes analysis within the General Meteorological Anal¬ 
ysis Package (GLMPAK) (des Jardins and Petersen 19S5). ranally, SGM analyzed 
temperature, relative humidity, and 1000 mb height fields, valid 6 h prior to each NOAA 
pass, were used as a first guess to compare with the other retrieval options. The retrieval 
process, with a regression first-guess and no surface data took about 15 min to run on 
a micro Yax 111 in the ID PA Lab. at NPS. The addition of surface data increased the 
processing time to about 30 min. The addition of SGM fields increased the processing 
time to nearly 1 h. 
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V. DATA ANALYSIS 


A. 0644 UTC 13 DECEMBER 19S8 
1. Overview 

The analysis of IOP-2 -starts with an evaluation of TOYS' data from the 
NOAA-11 pass just off the east coast of the U.S. at 0644 UTC 13 December 19S8. The 
initial IOP-2 cyclone was developing off the Carolina coast at this time. The NOAA-11 
pass (AYHRR, CH 4) (Fig. S) shows multi-layered cloudiness associated with the de¬ 
veloping low extending from the North Carolina coast eastward. Scattered low clouds 
are present northward along the coast from North Carolina to Cape Cod. with broken 
to overcast low clouds off the coast of Maine- Mostly high clouds extend from the Ohio 
Valley northward to Hudson Bay. These clouds are in advance of the upper level trough 
moving eastward into the Mid-Atlantic states. 

Seventeen ERICA rawimondcs and dropwindsondes were available for com¬ 
parison with collocated TOWS soundings (Tig. 9). Dropwindsondc and rawinsonde TOO 
mb winds were also plotted (Fig. 9a). A total of 600 soundings were retrieved 

from each NOAA-11 pass. For each case in this study, all dropwindsondes and 
rawinsondcs launched within 3-n of the IOYS soundings were used for comparison with 
the TOYS data. Dropwindsondes t .re named according to the aircraft type and time of 
launch. Those prefixed with an "A" were launched from L’SAF \YC-130's. while those 
prefixed with "N" were launched from NOAA P-3's. 'I he number codes represent the 
hour and nearest 10 min. they were launched. For example, A091 was launched at ap¬ 
proximately 0910 ETC. nearly 2 1 2 h after the NOAA-11 pass at 0644 UTC. TOYS 
soundings prefixed with an "11" are clear, while soundings prefixed with an ' F" arc cloudy 
(sec Chapter 4 for definitions of clear and cloudy TOYS soundings). Also, note that the 
TOYS soundings are not equally spaced, but represent the location of the warmest TOY 
in the 9-FOV box. So. some TOYS soundings overlap while others arc separated by a 
distance greater than 75 km. 

As discussed in the previous section, four different first-guess options were used 
to compare the TON’S retrieval data with the in situ data during IOP-2. They were: 
regression without surface temperature or dewpoint, regression with surface temperature 
and dewpoint, climatology with surface temperature and dewpoint, anti SGM temper¬ 
atures. relative humidity, and 10(><) mb neiglit fields with surface temperature and 





dewpoint (hereafter referred to as retrievals RN. RS, CS, and SS, respectively). The 
purpose of using all four methods is to understand the importance of the first-guess to 
the retrieval and to achieve the most accurate TOVS description of the ERICA cyclone 
environment. 

2. TOYS sounding comparisons with dropwindsondes and rawinsondes 

Three dropwindsondes, A091. AOS?, and A0S1, were launched in the cloudy 
baroclinic zone north and west of the initial IOP-2 cyclone. I he 0600 LTC rawinsondc 
for Bermuda (XKE) was available cast of the low center, while HAT and other coastal 
rawinsondes. available at the same time, were far to the west of the low. All four re¬ 
trieval methods for cloudy sounding. T454. aie compared to the XKF sounding. The 
regression with no surface (RN) and spectral with surface (SS) retrievals are presented 
in Tig. 10. while the regression with surface and climatology with surface retrievals ate 
shown in Fig. 11. 

The XKF temperature prolile shows a stable region from S5<)-7()<» mb that sep¬ 
arates weak easterly flow from the stronger westerlies aloft. The moisture trace indicates 
moist low level (100<>-850 mb) and upper level (5o()-?()() mb) conditions. '1 he RN re¬ 
trieval (Fig. 10a) failed to completely capture this stable structure, but the SS temper¬ 
ature retrieval (I'ig. 1 Ob) was very accurate. It resolves the stable temperaune structure 
from S5i»-70O mb and agrees well with the XKF temperature trace to 200 mb. However, 
th.e moisture profiles for both methods were too moist from 1 (mi<)- 7(K) mb and too dry 
above 700 mb. A comparison between F454 RS and C'S retrieval methods and XKF 
shows the RS and CS retrievals resolving the general temperature profile but unable to 
detect the low level inversion. All the moisture profiles indicated the low-level moisture 
but. in general, were poor. The XKF comparison clearly indicates the importance of the 
first guess to the retrieval. 

TOYS soundings were also compared with coastal rawinsondes in the mostly 
clear atmosphere w est of the l()P-2 cyclone. In general, the TOYS SS retrieval method 
produced accurate temperature profiles compared with the rawinsondes. as found with 
XKF. I or example. HAT was compared with the collocated TOYS sounding. F424. in 
Fig. 12. The difference between the F424 RN and SS retrievals with HAT is striking. 
The RN retrieval misses the cold low-level layer and associated inversion from 900-750 
mb. In contrast, the SS sounding does resolve the general structure of the low- 
troposphere. However, both retrievals miss the sharp tropopamc at 35o mb in the HA T 
sounding. Although the RS and. CS retrieval methods for F424 (not shown' were not 
quite as accurate as the SS retriesal. they did accurate!) depict the low level tempo; ature 
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structure because of the inclusion of surface temperatures in the first guess, flic 
moisture profile from 100u-500 mb was fairly well represented by all four retrieval 
methods. However, the shallow dry layer near 850 mb was not resolved. It is unlikely 
that this dry layer would be captured by the coarser vertical resolution of TOYS. 

A comparison between other coastal rawinsondes and collocated TOYS 
soundings produced similar results to the HAT case. The TOYS soundings w : th the SS 
retrieval method captured most of the thermal structure measured by the rawinsondes. 
The RS and CS methods, although not as accurate as the SS method, were able to ac¬ 
curately capture the colder air from 1000-700 mb. The RN method was not successful 
in resolving the low level inversion. The collocated TOYS soundings using all four 
methods gave successful retrievals compared to rawinsondes located further inland (PIT 
and FNT). These soundings had a fairly constant lapse rate and were easier for TOYS 
to retrieve. However, the SS retrievals remained the most accurate for temperature and 
moisture. 

Over the cloudy region associated with the initial IOP-2 cyclone, the TOYS 
sounding comparisons were not as successful. Dropwindsonde A0S1. at 35°N 74°\Y, 
was launched in the western edge of the baroclinic zone associated with the developing 
cyclone. It shows stable frontal structure in the low levels and wind shear similar to tire 
XKF case (Tig. 13). The A0S1 sounding is saturated to 500 mb and suggests that rain 
is probably falling in the vicinity. TOYS sounding L426 was used for comparison with 
A0S1. The RN retrieval (not shown) failed to capture any of the temperature structure 
of A0S1. Instead, it showed a smooth temperature curve throughout the troposphere, 
3°C too warm at 1000 mb and 3°C too cold at 500 mb. The RS 13426 retrieval (I'ig. 13a) 
shows a more accurate profile from 1000-800 mb than the RN method, but fails to show 
the inversion from 850-700 mb. The CS method (not shown) produced a profile identical 
to the RS method. The SS retrieval (fig. 13b) does suggest the inversion structure from 
850-700 mb. Although it is too cold above 850 mb. it represents the weak temperature 
inversion better than the other options. All four retrieval methods derived useful 
moisture profiles indicating the saturated state of the A0S1 profile. One other compar¬ 
ison (not shown) was made between A0S1 and a collocated TOYS retrieval. 1 427, a 
cloudy sounding just southeast of A0S1 (sec fig. 9), was quite similar to A081. The SS 
retrieval resolved the weak inversion, while the other retrieval methods produced smooth 
temperature profiles. 

The remaining overwater dropwindsondcs, A09I and A0S3 (not shown), were 
located closer to the surface low. Both were saturated up to 500 mb. indicating precip- 
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itation in the region. They also showed a stronger stable region in the low troposphere, 
similar to the XKT sounding. The collocated TOYS soundings using all four retriesal 
methods were too warm at 1000 mb and too cold at 500 mb. However, the SS retrievals 
were still the most accurate. A comparison between all four TOVS retrieval methods 
with the dropwindsondes and rawinsondes in the region indicated that only the SS re¬ 
trievals svere successful in capturing some of the structure of the cloudy atmosphere 
shown by the in situ data. The RN retrievals gave smooth profiles, attempting to rep¬ 
resent the average temperature in the troposphere and lower stratosphere. The RS and 
CS retrievals produced more accurate temperature profiles from 1000-700 mb. but still 
failed to capture the thermal structure in the lower troposphere. However, even though 
the SS retrievals were the most accurate, only the TOVS CSS; retrieval collocated with 
XKT depicted all of the cloudy temperature structure shown by the rawinsonde. 

The TOYS overland retrievals were much better than the overwutcr retrievals 
for this pass. However, the land area was dominated by high pressure and mostly clear 
skies with a constant lapse rate for most soundings. T hus, the temperature profile was 
considerably easier to recover from the satellite data. Over water, where cyclogcnesis 
was occurring, the TON’S soundings were less accurate compared with the 
dropwindsonde measurements. Rain appeared to attenuate most of the TON’S 
soundings in the region, making them too cold at 500 mb. As a result, the inversion 
structure of the dropwindsondes was not totally resolved in the TON'S soundings. The 
use of SGM analyses for the first guess, even though 6 h old. improved the TON’S tem¬ 
perature structure considerably. The surface data were critical to TOYS resolving the 
boundary layer structure, particularly for the coastal rawinsondes. 

3. Cross section comparisons 

Another useful method in diagnosing cyclogcnesis is constructing potential 
temperature cross sections through interesting regions of the cyclone. Significant com¬ 
ponents of the cross sections include steeply-sloped or packed isentropes and wind speed 
and direction changes to identify the location and strength of fronts. The ION’S data 
do not include actual wind information, but location and intensity of fronts and the 
tropopausc level can be determined by the slope and packing of the isentropes. 

f or this first case, a cross section was constructed using HAT, A081, A0S3, 
A091, and XKT through the northern quadrant of the initial IOP-2 cyclone (fig. 9a). 
Tig. 14 compares this cross section with a collocated oserwater ION'S cross section us¬ 
ing the SS retrievals. The observed cross section located along the warm front, and the 
packed isentropes in the region between AOS3 and XKT indicate the strongest part of 







the front. Tire TON'S cross section (Fig. 14b) resolves some of the frontal structure. 
Also, note the steeply sloped isentropes between L424 and F427, showing the cold air 
along the coast similar to the isentropes between HAT and A0S1. Even though the 
TOYS cross section isentropes are not as tightly packed as observed farther cast (be¬ 
tween F430 and F454), they still indicate the strongest area of the baroclinic zone fairly 
well. The tight gradient of the isentropes from 300-200 mb also resolves the location of 
the tropopausc quite well. The TOYS tropopausc corresponds to the levels of maximum 
winds at HAT and XKF. 

4. Static stability comparisons with SGM fields 

The the static stability parameter, , described by Chaifant (19S9) in a previ¬ 
ous study of IOP-2, was evaluated for the TOYS data set. Low static stability provides 
a more favorable environment for upward vertical motion and cyclogenesis. During the 
initial stages of IOP-2. at least two cyclones were present. The first cyclone, analyzed 
by Chaifant to be east of Florida at this time, was forecast to rapidly develop by the 
NMC and NORA PS models (Miller 1990). But a second low developed to the northwest 
of the initial low. It was this system that rapidly deepened from OO00-120O E TC 14 
December I9SS. 

Chaifant found a minimum in 500-1000 mb. in the NMC SGM analvscs 
.oir the Carolina coast for 1200 ETC 13 December (Fig. 15b). The second cyclone de¬ 
veloped in this region during the next 12 hours. The dropwindsondes also confirmed this 
lower stability near the coast. The AOS I profile was more unstable than the coastal 
rawinsondcs to the west and the dropwindsondes to the cast. Stability from the TON'S 
0644 ETC data was also evaluated to determine if the areas of low static stability de¬ 
scribed by the SGM analyses were also present in the TON'S data. The TOYS stability 
analysis at 0644 E TC (Fig. I5at shows lower static stability olf the North Carolina 
coast, in the region where the eventual rapid deepener began to develop at 1 S 00 ETC. 
This is in agreement with the SGM analysis that shows the lower static stability over a 
much broader region. Note that operational TON'S soundings would be included in the 
SGM analyses via the data assimilation cycle. So it is not surprising that the two ana¬ 
lyses arc similar. Also, note that the TON'S stability analysis has less meaning on the 
eastern edge, as the TOYS data did not extend that far cast (Fig. 9 1 . 

Although both the TOYS and SGM regions of low static stability are partially 
confirmed by the dropwindsondes and indirectly by subsequent surface cydogenesis. the 
spatial distribution and absolute magnitude is somewhat uncertain. This is partially due 
to the tendency of the cloudy TON'S soundings to be too warm at the surface and too 






cold aloft, effectively lowering the static stability. While spatial smoothing removes 
some of this error, a future study to determine the error distribution and its relationship 
to cloudy retrievals would be illuminating. 

5. Height and temperature comparisons with SGM fields 

One final comparison was made between TOYS and SGM gcopotential heights 
and temperatures. Of interest was the quality of the TOYS (plus surface data) 
gcopotential height and temperature analyses for this initial cyclogenesis period. All 
mandatory levels from 1000 mb to 500 mb were evaluated. Only 7o0 mb level compar¬ 
isons will be shown, as this level proved to be the most accurate. However, the inclusion 
of a 1000 mb height analysis in the first guess is critical to the retrieval of accurate 
heights at all mandatory levels. Otherwise, less accurate climatological heights are used, 
which produce poor heights at all levels. 

f ig. 16 compares the TON'S Too mb height and temperature analyses at 06-14 
UTC with the SGM analyses at 1200 l.'TC. The TOYS 700 mb analyses (Pig. 16a/ show 
the weak trough associated with the first cyclone off the North Carolina coast. The 
height analysis also shows general ridging over New Iingland. The TOYS analyses are 
in general agreement with the SGM analyses, 6 h later (Tig. 16b). 

In conclusion, the overall effectiveness of TON'S data at 0644 LTC was very 
encouraging, especially with the inclusion of SGM fields in the first guess anJ the use 
of a surface temperature analysis. The TOYS cross section using the SS retrieval method 
resolved the warm frontal region associated with the first I OP-2 low. The TOYS sta¬ 
bility analysis isolated the lower static stability over the coastal waters, which was con¬ 
firmed by the dropwindsondc data. The TON’S analyses did resolve the 7O0 mb trough 
associated with the developing cyclone. Because the SS retrieval method was so superior 
to the other methods, only the SS retrieval will be shown for comparison with the other 
data sources for the remainder of this thesis. 

B. 1820 UTC 13 DECEMBER 1988 

Twelve hours later a second NOAA pass covered the western North Atlantic at 1S20 
UTC. By this time, the second surface low had developed cast of N. Carolina, near 
37 C N 72°\V, in the relatively unstable air northwest of the initial low. The AY11RR 
imagery (Fig. 17) shows the cloudiness with both cyclones and strong convection around 
and south of the new low. The imagery reveals the rapid cloud growth of the second low 
in response to the movement of an upper level short wave trough over the coastal area. 



Twenty-eight rawinsonde and dropwindsondes were collocated with the TOYS 
soundings from this NOAA-11 pass (fig. IS). Dropwindsondc and rawinsonde 7()0 nth 
winds were also plotted (Fig. ISa). Five dropwindsondes (N14U. X’142. \144. N' 150 . and 
"N152) were launched at approximately 1500 FTC. Their locations arc indicated on an 
AVURR CII 1 subscene (Fig. 19). The evaluation of the TOYS retrievals at 1820 FTC 
starts with comparisons in this convective region near the old frontal zone, but ahead 
of the intense convection of the developing 10P-2 storm (see Fig. 19). The satellite im¬ 
agery shows partly cloudy conditions in the vicinity of X150 and X144. Their locations, 
along with collocated TON’S soundings. E316 and F3IS. arc displayed in f ig. 20. The 
only significant variation in the vertical structure shown by the dropwindsondes is a 
weakly stable layer from 850-7OO mb by XI50 and from about 6 S 0 - 6 OO mb by X144. 
The temperature profile for the TON’S sounding, K316. was very accurate compared with 
XI50 (Fig. 20a). The ION’S profile was within FC of the dropwindsondc from 1000-50'> 
mb. except for a slightly larger departure near 900 mb. In the comparison between X144 
and FA 18 (Fig. 20bi. F 3! S was not quite as successful. It was about 5‘C too warm at 
1000 mb, but within 2~C at 850 . Too. and 5uo mb. The weak inversion at 6S0-600 mb 
is p.ot resolved by the TON’S sounding. And the large temperature discrepancy at 1000 
mb may be due to inaccuracies in the surface temperature analysis used in the retrieval 
or because of the time difference between X144 and F 518 (at least 3 hi. 

Two other comparisons were made further north along the frontal zone in cloudier 
areu.v In these areas, the dropwindsondc temperature profiles, X142 and XT40. showed 
a stronger inversion than the two previous profiles. X142 and X14o were compared with 
TON’S profiles. F3oo and L2S2. respectively (Fig. 21). In both cases , the TON’S profiles 
were too warm from 1 000 - TOO mb and too cold above 700 nib. Both dropwindsondes 
show a stable layer near Too mb that is too shallow to be captured by the I ON’S re¬ 
trievals. However, the F2S2 temperature profile (Tig. 21b) shows some indication of 
attempting to represent the stable layer, but only weakly. 

Other comparisons between overwater dropwindsondes and TOYS soundings 
showed similar results. The TON’S soundings were fairly accurate, but were the least 
accurate where shallow inversions existed between mandatory levels. The TON’S 
soundings were more accurate when the inversion layer was deeper and best described 
the soundings with a constant lapse rate. The TOYS moisture profiles did not compare 
favorably with the dropwindsondes as. in general, they were too moist. 

The comparisons made between TON’S soundings and coastal rawinsondcs showed 
less success than earlier. For example, the HAT and NY A L. temperature profiles were 



much different than their collocated TOYS profiles, F204 and F257 (fig. 22), even 
though the atmosphere was mostly clear. The tropopausc lowered to about 400 mb at 
both locations due to the strong upper level trough over the region (Refer to previous 
figure of HAT at 0644 UTC). The TOYS retrievals analyzed the tropopausc at 300 mb, 
and the resulting TON'S temperature profiles were 4-S c C too warm from about 700-300 
mb for both cases. The lower level temperature structure was fairly accurate though. 

The coastal rawinsondcs and TOYS profiles further north compared more favorably. 
For example, the ACY and Cl III temperature profiles were very well represented by the 
collocated TOYS soundings, E221 and II205 (not shown). The F22I temperature profile 
was within 2°C accuracy of the ACY profile front 1000 mb to about 550 mb, except for 
a very shallow inversion layer at 700 mb shown by tltc ACY profile. The tropopausc 
level of ACY was at about 350 mb. while F221 indicated a 300 mb tropopausc, mini¬ 
mizing the upper troposphere errors. The tropopausc at CH11 was at 250 mb, and the 
E205 tropopausc matched it almost perfectly. The result was an almost perfect match 
between the two soundings ( < l'C error), except for the layer from 1000-S50 mb. where 
the TOYS profile was warmer. 

A cross section in the region of the dropwindsondcs launched at 1500 UTC (Fig. 
23a). shows the frontal structure between the two surface centers. The collocated TOYS 
cross section (Fig. 23b) reproduced this thermal pattern. The slope of the isentropes in 
the two cross sections is nearly identical. Also, the packing of the TON'S low 
troposphere isentropes indicates the location of the front in the cross section 600-$<K) 
mb layer. The tropopausc also is easily identified on the TOYS cross section. 

Another TOYS cross section was constructed across the large-scale barodinic zone 
associated with the developing low, from F.221 (4<f N 74 C NY), northwest of the cyclone 
along the New Jersey coast, southeast to F357 (33 c \ 67 C NY) in the warm sector. The 
cross section (Fig. 24) depicts the thermal structure of the developing cyclone. It shows 
steeply sloped isentropes just cast of the low and more tightly packed isentropes associ¬ 
ated with the large-scale barodinic zone and the tropopausc aloft. The weak vertical 
stability is also evident in the lowest 200 mb southeast of sounding l;200. 

TOYS static stability from the 1S20 UTC 13 December pass (I ig. 25a) was com¬ 
pared with the SGM stability at 0000 UTC, 14 December (Fig. 25bI. 1 he TOYS stability 
analysis continues to indicate a stability minimum near the center of the developing 
cyclone. The SGM stability field shows a larger region of weak stability. Unfortunately, 
the dropwindsondcs arc cast of the TOYS minimum and cannot confirm this feature in 



cither analysis. However, this region of w eak stability docs correspond to the region of 
subsequent vigorous surface development, which is dynamically consistent. 

The TOVS 700 mb temperature and height analyses for IS 10 13 December arc now 
compared to the SGM 0000 L’TC 14 December analyses (Tig. 26). The TOYS height 
analysis successfully captured the developing trough cast of North Carolina associated 
with the rapid cyclogcncsis event. Also evident is the warm advcction region located in 
the vicinity of the growing comma cloud (Tig. 19). The trough location is consistent 
with the 0000 L’TC SGM analysis which shows eastward movement and continued am¬ 
plification. The TOYS isotherm pattern captures the 700 mb frontal zone off the Last 
Coast, but fails to show the characteristic "S" shape isotherm pattern of a developing 
disturbance. 

In conclusion, the TOYS evaluation at 1S20 L'TC showed that the TOYS data were 
very useful in diagnosing the details of the IOP-2 cyclone. The collocated TOYS 
soundings continued to compare favorably with the dropwiudsondcs. However, the 
TOYS soundings were the least accurate in regions of strong inversions or thick 
cloudiness, analyzing the surface too warm and the mid levels too cold. The sounding 
temperatures were the most accurate at about 7<>0 mb. The only significant limitation 
of the TOYS soundings was their inability to capture the lowering tropopausc in the 
region of the upper level trough along the coast. The low level temperature profiles were 
successful in nearly all of the comparisons. 

C. 0637 L'TC 14 DECEMBER 1988 

Fxplosive development of tire extratropical cyclone was now in progress at 0637 
L'TC, the tunc of the next NOAA-11 pass. The storm’s central pressure dropped about 
30 mb from OODD-1200 LIT, as it moved eastward (Chalfant 1989). TRICA aircraft 
launched three dropwiudsondcs around its center near the time of this pass. Fig. 2 m 
shows the dropwindsondcs and rawinsondes launched within 3 h of the NOAA-11 pass 
at 0637 L’TC, along with 700 mb winds. Fig. 27b shows the collocated 1’OYS soundings. 
One of the dropwiudsondcs, N023 (not shown), failed, only sampling the atmospheric 
layer from about 700-500 mb. The other two dropwindsondcs. N020 and N070, were 
used for comparison with collocated TOYS soundings (Fig. 28). 

The satellite imagery over the region (not shown) shows thick cloudiness over the 
area. However, the TOYS sounding. FI82. collocated with N020 (Fig. 28a), produced 
a fairly accurate temperature profile, compared with the dropwinJsonde. It was ap¬ 
proximately 3 ( too warm at the surface and missed a shallow inversion layer between 



700-500 mb, but otherwise produced a useful temperature profile. The TOYS sounding, 
F204, collocated with N070 (Tig. 28b). produced a fairly accurate profile, but it was 
3-5°C too cold throughout the lower and mid troposphere. This cloudy TOYS sounding 
is cold likely due to the effects of precipitation attenuation of the MSU channels. 

Other comparisons between rawinsondcs and TOYS soundings for this time period 
(not shown) produced similar results to the previous time periods. The TOYS temper¬ 
ature profiles were generally within 2-3’C of the rawinsonde profiles at all the mandatory 
levels in the troposphere. The largest errors were at 1000 mb and the tropopause. Un¬ 
fortunately, no conclusions can be made about their accuracy near an explosively de¬ 
veloping cyclone, because only two comparisons could be made. However, additional 
insight can be gained by comparing TOYS cross sections and objective analyses over a 
larger area. 

A cross section of the 0637 UTC TOYS soundings (Fig. 29) was constructed through 
the warm front of the rapidly developing cyclone, from F143 (41 5 N 67°\Y) north of the 
low center near Cape Cod, southeast to F321 (30°N 63°NY) in the warm sector. The 
cross section shows steeply sloped isentropes ahead of the warm front, but misses the 
location of surface front near F.244, according to Chaifant's (1989) analysis. However, 
the cross section shows tighter isentropes associated with the warm front aloft. Also 
note the sloping tropopause and location of the jet stream clearly defined by the TON’S 
cross section. 

The minimum stability from the TON’S analysis (Fig. 30a) was again located near 
the center of the cyclone. The I ON'S analysis was compared with the SOM analyses at 
0000 UTC, 14 December (Tig. 25b) and 1200 UTC, 14 December (Tig. 30b). Notice the 
similarity of the TON’S analysis to the SGM analysis at 0000 UTC. However, the TON’S 
analysis places the lower stability area northeast of the SGM analysis, and shows higher 
stability southwest, where the SGM minimum stability area is located. Recall that these 
two analyses arc not independent because TOYS data are used in the SGM analysis. 
At 1200 UTC, the SGM analysis concentrates the minimum stability area into a more 
confined location, in the warm sector south of the position of the surface low at 1200 
UTC 14 December 1988. '1 he TON’S 0637 UTC minimum stability area is in the ap¬ 
proximate location of the surface low at 1200 UTC. Some of these differences between 
the SGM and TON’S static stability analyses may be due to the retrieval problems in 
cloudy areas as well as differences in time of analysis. This point should be explored 
more thoroughly in a future study. 









The TOYS 700 mb height and temperature analyses (Fig. 31a) at 0637 FTC 14 De¬ 
cember indicates continued development of the 700 mb trough. The 0000 FTC SGM 
analyses (Fig. 26b) shows a strong trough developing just cast of the mid-Atlantic coast. 
By 1200 FTC (Fig. 31b). it has developed into a cutolflow near 35°N 65TV. The TOYS 
analysis deepens the trough more than the 6-h earlier 0000 FTC SCi.M analysis. The 
TOYS temperature analysis looks fairly reasonable compared with the SGM analyses, 
but may be too cold north of the low. Still, the TON’S analysis shows strong warm 
advcction occurring ahead of the low. consistent with vigorous cyclogencsis. 

D. 1810 FTC 14 DECEMBER 1988 

The analysis of IOP-2 ends with an evaluation of TOYS data at 1S10 FTC 14 De¬ 
cember 1988. The NOAA-U AVHRR. imagery (trot shown) portrays a well-wrapped 
cyclone at the far eastern edge of the pass near 3S C N 63°NV. Many dropwindsondes were 
launched within 3 h of the NOAA-11 pass around the storm center. Their locations and 
7t.)0 mb wind observations arc depicted in Fig. 32a. Note how well the dropwindsondc 
winds delineate the center of the 7(H) mb trough. Unfortunately, the region of TON’S 
soundings (Fig. 32b) onlv -•:* nds to the cyclone center. The intense frontal /tone cast 
and north of the low i c v .* covered by this pass. 

Dropwindsondes, M71 and N’203, were compared with TOYS soundings, L205 and 
E225. respectively (Fig. 33). In both cases, the TOYS soundings represented the 
dropwindsor.de temperature profiles reasonably well. The 'I ON’S soundings were colder 
than the dropwindsondes likely due to precipitation in the region. The SGM first-guess 
fields in the TON'S retrievals were extremely important for retrieving accurate TON'S 
profiles, as seen at all the previous times. 

A cross section (Tig. 34) from F149 (44°N 63TY) north of the storm near the Nova 
Scotia coast, through the low to F263 (37°N 62 3 NY) along the cold front, showed the 
location of the surface trough (denoted by the slope of the isentropes near FI87) and a 
baroclinic zone north of the low. It also depicted the location of the tropopausc quite 
well. In summary, the TON’S analyses at 1810 FTC were still very useful in describing 
the cyclone environment. TOYS gcopotential height and temperature analyses were not 
presented for this case because the 'ION’S data covered only the western half of the 
cyclone circulation. 



VI. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

This thesis examined the effectiveness of TOVS in describing the environment of the 
IOP-2 cyclogenesis event. Numerous dropwindsondes and IOVS comparisons were 
made to assess the ability of TOYS to retrieve temperature and moisture profiles. Em¬ 
phasis was placed on the synoptic evaluation of the TOYS soundings and associated 
horizontal fields rather than computing error statistics. TOYS soundings were found to 
be effective in describing the synoptic-scale distributions of static stability and upper 
level height and temperature fields associated with a developing cyclone. TOYS poten¬ 
tial temperature cross section^ also successfully described the larger-scalc structure of the 
cold frontal or warm frontal region associated with the iOP-2 cyclone, as well as the 
location of the tropopause. 

However, a limitation in the overall effectiveness and use of TON'S soundings is that 
TOYS retrievals require an accurate first guess to derive accurate profiles. In the si¬ 
multaneous retrieval method, surface data in the first guess proved to be crucial in the 
accurate retrieval of temperature Irani the surface to 700 nib. And. the addition of 
model gridded temperatures was important for retrieving vertical structure in the lower 
troposphere. The further addition of loop mb heights in the first-guess was crucial for 
obtaining a physically reasonable height field from the TON'S data, finally, an impor¬ 
tant observation of this study was that the TOYS overland soundings generally were 
more accurate than the overwater soundings. A primary reason for this was that most 
of the overwater soundings were in areas associated with cyclogcncsis during IOP-2. 
The least accurate soundings were those affected by clouds and or precipitation. 

B. RECOMMENDATIONS 

Further study of TOYS data should include the use of forecast fields in the first 
guess instead of 6-h old observations. Now that an interface is available to input model 
gridded data in the IDEA Lab., comparisons can easily be made between other gridded 
data sets, such as the NORAPS forecast fields or NGM forecasts, to compare the ac¬ 
curacy of the retrievals with different first-guess methods. The use of gridded data 
greatly enhanced the TON'S retrievals in this study, but the inclusion of surface data was 
also crucial. Unfortunately, the additional surface data included in the TON'S retrievals 
for thh stud;, will not be available operationally. 





Another option to test is the use of the National Lnvironmcntal Satellite Data and 
Information Service (N'LSDIS) operational gridded sea-surface temperatures derived 
from multi-channel NOAA AY1IRR. This option was already explored by NOARI. 
scientists (Richardson and Haugen 19S9). Sea-surface temperatures were input into the 
simultaneous retrieval as surface skin temperatures (versus first-guess surface temper¬ 
atures) during the Genesis of Atlantic Lows Lxperimcnt (GALL). Their inclusion 
produced more accurate overwater surface temperatures compared to a NORAPS sur¬ 
face analysis. Additionally, the sea-surface temperatures improved low level cold frontal 
structure on cross sections. However, very little difference (with and without sea-surface 
temperatures) was observed in the temperature profile above 1000 mb. The study con¬ 
cluded that the use of sea-surface temperatures in the retrieval needs to be studied fur¬ 
ther. 

Several options described in the original Lxport package were not available for in¬ 
corporation into this thesis. They include: the filtering of redundant or inaccurate 
TON'S retrievals before they are used in objective analyses; the derivation of geostrophic 
winds based on the temperature and geopotential height retrievals; and the availability 
of higher resolution geography in the determination of surface elevation. As discussed 
in the previous chapter, a study examining the error distribution of TON'S soundings in 
cloud vs. no-cloud retrievals would be helpful in making better use of TON'S soundings 
in weather forecasting. The filtering option, mentioned above, should be incorporated 
into such a study to determine its significance. And finally. TOYS MSI.' channel images 
might provide some useful information in determining the location of low pressure cen¬ 
ters and frontal regions and should also be evaluated in future TON’S studies. In con¬ 
clusion. TON'S data provide a useful supplement to other data sources and should be 
incorporated into future cyclogencsis studies. 
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Tig. 3. TOYS and radiosonde composites. 1 he rms diH'crcncc between North 
American interactively derived IIROS-N soundings and radiosondes. 22 March-19 
April 1979 (I rom Smith ct al. 1979). 
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I" jo. 5. Comparisons of 1 ()\ S and UCMWU 1000-51*0 mb thickness 
analvscs. 1 he lOYS analyses (I ig. 5a) shows a higher amplitude ridge between 
Uneland and Spain and a deeper trough over southern Europe than the UCMWU 
analysis (Fig. 5b) (From Smith ct ah 1983). 
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Fig. 8. NOAA-11 AVHRR CH 4 imagery, 0044 L'TC 13 December 19SS. 
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J ig. 9. Rawinsonde and dropwindsondc locations at <»<>»K> L 1C and 
collocated IOYS soundings at 0644 C[C 13 December 19SS. lie. 9a shows 
rawinsonde and dropwindsonde locations, f ie. 9b shows JOYS sounding locations. 
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l ie. 13. Dropwindsoiulc AUS1 with collocated K)VS sounding 1 ig. 13a 

compares AbSI with a RS lOWS iciiicval. l ie. 13b compares AOS1 with a SS re¬ 
trieval. 
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I ip. 14. Cross section through rawinsondes and overwater dropwindsondes, and 
collocated TOVS cross section. 1 ip. 14a shows the cross section with the in situ 
data. l ie. 14b shows the cross section with the collocated I OVS soundines. 


































l ie. 15. Comparison between I OX'S static stability analysis and the NMC oper¬ 
ational analysis I'rom the SGM (> It later. I ie. 15a shows the I U\ S static stability 
analysis at 0644 l."IC 13 December I9SS. 1 je. 15b shows the NMC analysis at l2oo 
L IC 13 December 19NS. 







I'ig. 16. Comparison between TOYS 7uo mb height and temperature anahscs and 
the N.MC operational analyses from the SOM 6 h later. Tig. 16a shows the lOYS 
analyses at 0644 UTC 13 December l'OS. l ie. 16b shows the NMC analyses at 
I2<)() L"IC 13 December 1988. 
























Fig. 18. Rawinsondc and dropw'inJsomJc locations at 15<M)-2)00 I K' and 
collocated lOVS soundings .it 1S2(» L I (J 13 December 1988. I ig. 18a shows 
rawinsondc and dropwindsondc locations. l ie. 18b shows IONS sounding lo¬ 
cations. 
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lie. 20. Dropuindsondcs N15o and NI4-J with collocated ION'S profiles 1016 
and l:31S. I ie. 2Ua shows the diflciciKe between N150 and 1:316. l ie. 2ob shows 
N144 and 11318. 
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I'ig. 21. Dropwindsondcs N142 and NMD with collocated 1 OVS profiles 1.300 
and L2S2. l ie. 21a shows the diHercncc between N142 and U300. I ig. 21b shows 
N140 and 12282. 
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Fig. 23. Dropwindsondc and IOYS cion sections at 1820 1 1C. l ie. 23a shows 
a cross section through the drnpwindsondcs at 1500 l.' K'. F ig. 23b shows a cross 
section through the collocated iOYS soundings. 































































Fie. 26. Same as f ig. 16 except for M>YS data at 1820 l'IC 13 Deember 1988 
and NMC data at 0000 UTC 14 December 1988. I ig. 26a shows the lOYS 700 
mb analyses. Fig. 26b shows the SC1M TOO mb analyses. 










l ie. 27. Rawinsonde and dropuindsonde locations at OhiO-OWO l J C' and 
collocated TONS soundings at 0637 I K' 14 December 19XS. lie. 27a shows 
rawinsonde and dropwindsonde locations. Fig. 27h shows TOYS sounding lo¬ 
cations. 

















1'ig. 28. Dropwindsondes N020 and No7<> with collocated IOVS profiles I' 1S2 
and 1’2‘M. I ig. 2Sa shows the difference between NU2U and 1:182. l ig. 2Sb shows 
N070 and F20-4. 








. 29. lOVS cross section over developing cyclone at 063™ I K' 14 December 
S. Cross section extends from I 143 (41 N 67 5 \V) to 1:321 (30'N 63'\V). 





















































lie. 32. Rawinsonde and dropwindsondc locations at I»*(> l.'IC and 

collocated ION'S soundings at IS|o I. 10 14 December 19SS. I ig. 32a shows 
rawinsonde and dropwindsondc locations. I ig. 32b shows ION’S sounding lo¬ 


cations. 






















lie. 33. Dropwirulsoiules N17I and N2o3 with collocated ION'S profiles I;205 
and 1:225. i ic. 33a shows the difference between N17J and 132o5. lig. 33h shows 
N203 and 1:225. 
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